Estimation and Prediction for Spatial Generalized Linear Mixed Models Using High Order Laplace Approximation
نویسندگان
چکیده
Estimation and prediction in generalized linear mixed models are often hampered by intractable high dimensional integrals. This paper provides a framework to solve this intractability, using asymptotic expansions when the number of random effects is large. To that end, we first derive a modified Laplace approximation when the number of random effects is increasing at a lower rate than the sample size. Secondly, we propose an approximate likelihood method based on the asymptotic expansion of the log-likelihood using the modified Laplace approximation which is maximized using a quasi Newton algorithm. Finally, we define the second order plug-in predictive density based on a similar expansion to the plug-in predictive density and show that it is a normal density. Our simulations show that in comparison to other approximations, our method has better performance. Our methods are readily applied to non-Gaussian spatial data and as an example, the analysis of the rhizoctonia root rot data is presented.
منابع مشابه
Parameter Estimation in Spatial Generalized Linear Mixed Models with Skew Gaussian Random Effects using Laplace Approximation
Spatial generalized linear mixed models are used commonly for modelling non-Gaussian discrete spatial responses. We present an algorithm for parameter estimation of the models using Laplace approximation of likelihood function. In these models, the spatial correlation structure of data is carried out by random effects or latent variables. In most spatial analysis, it is assumed that rando...
متن کاملBayesian Inference for Spatial Beta Generalized Linear Mixed Models
In some applications, the response variable assumes values in the unit interval. The standard linear regression model is not appropriate for modelling this type of data because the normality assumption is not met. Alternatively, the beta regression model has been introduced to analyze such observations. A beta distribution represents a flexible density family on (0, 1) interval that covers symm...
متن کاملPractical likelihood analysis for Spatial Generalized Linear Mixed Models
We propose a standard approach to make inference for spatial generalized linear mixed models using Laplace approximation. Based on analysis of two datasets previous analysed in literature, we compare our approach with different approaches. The first the rhizoctonia root rot dataset is an example of Binomial SGLMM and the second rongelap dataset is an example of Poisson or Negative Binomial SGLM...
متن کاملComparison of Different Estimation Methods for Linear Mixed Models and Generalized Linear Mixed Models
Linear mixed models (LMM) and generalized linear mixed models (GLMM) are widely used in regression analyses. With the variance structure dependent on the random effects with their variance components, the parameter estimation of LMMs is more complicated than linear models (LM). Generally, we use maximum likelihood estimation (MLE) together with some procedure such as derivative free optimizatio...
متن کاملLaplace approximation for measurement error correction in generalized linear (mixed) models
The Laplace approximation is applied to measurement error models when the likelihood function involves high dimensional integrals. The cases considered are generalized linear models with multiple covariates measured with error and generalized linear mixed models with measurement error in the covariates. The asymptotic order of the approximation and the asymptotic properties of the Laplace-based...
متن کامل